Phase synchronization in cerebral hemodynamics 
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A healthy human brain is perfused with blood flowing laminarly through cerebral vessels, pro- 
viding brain tissue with substrates such as oxygen and glucose. Under normal conditions, cerebral 
blood flow is controlled by autoregulation as well as metabolic, chemical and neurogenic regulation. 
Physiological complexity of these mechanisms invariably leads to a question as to what are the re- 
lations between the statistical properties of arterial and intracranial pressure fluctuations. To shed 
new light on cerebral hemodynamics, we employ a complex continuous wavelet transform to deter- 
mine the instantaneous phase difference between the arterial blood pressure (ABP) and intracranial 
pressure (ICP) in patients with traumatic brain injuries or spontaneous cerebral hemorrhage. For 
patients with mild to moderate injury, the phase difference slowly evolves in time. However, severe 
neurological injury with elevated ICP are herein associated with synchronization of arterial and in- 
tracranial pressure. We use Shannon entropy to quantify the stability of ABP-ICP phase difference 
and discuss the clinical applicability of such measure to assessment of cerebrovascular reactivity and 
autoregulation integrity. 

PACS numbers: 87.80.Tq, 87.19.Uv, 87.10. +e 



Synchronization between different physiological sys- 
tems or subsystems has long been recognized as a ubiq- 
uitous dynamical effect [l|. Examples of such synchro- 
nization are as diverse as circadian rhythm jjj, |2|], cor- 
relation of respiration with mechanical ventilation [3| or 
locomotory rhythm [J], coordinated motion [l|, animal 
gait |ij], or synchronization of oscillations of human in- 
sulin secretion with glucose infusion J6J. These phenom- 
ena are essentially confined to nearly periodic rhythms. 
However, the development of novel concepts [3, lsj| has 
paved the way for the application of synchrony analysis 
to inherently nonstationary and noisy signals; time series 
that are characteristic of cardiology and encephalography 
lllllil (see also references therein). 

A healthy human brain is perfused with blood flow- 
ing laminarly through cerebral vessels, providing brain 
tissue with substrates such as oxygen and glucose. Cere- 
bral blood flow (CBF) is relatively stable, with typical 
values between 45 and 65 ml/lOOg of brain tissue per 
second, despite variations in systemic pressure as large 
as 100 Torr. This phenomenon, known as cerebral au- 
toregulation (CA) [llj, is mainly associated with changes 
in cerebrovascular resistance of small precapillary brain 
arteries. CBF is also affected by metabolic, chemical and 
neurogenic regulation. Strong susceptibility of brain tis- 



sue to even short periods of ischemia underlies the physio- 
logical significance of these intricate control mechanisms. 
In the phenomenological description of cerebral hemo- 
dynamics, fluctuations of arterial blood pressure (ABP), 
due to pressure reactivity of cerebral vessels, lead to fluc- 
tuations of intracranial pressure (ICP). The goal of this 
paper is to quantitatively analyze the interplay of ABP 
and ICP from the viewpoint of synchronization. In par- 
ticular, we investigate time evolution of instantaneous 
phase difference between ABP and ICP time series in 
patients with traumatic brain injuries or spontaneous 
cerebral hemorrhage. We examine the extent to which 
the relative phase dynamics reflects pathological condi- 
tions. We adopt the mathematical framework of synchro- 
nization since it is particularly well-suited to the analy- 
sis of non-stationary bivariate time series. Interestingly 
enough, we have not been able to find previous applica- 
tion of this a ppr oach to cerebral hemodynamics, see, for 
example 
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Let us consider two signals Si(t), s 2 {t) and their corre- 
sponding instantaneous phases 4>\{t) and <f>2(t). Phase 
synchronization takes place when n<f>i{t) — m<f>i(t) = 
const where n, m are integers indicating the ratios of 
possible frequency locking. Herein we consider only 
the simplest case n = m = 1. Furthermore, as with 



most biological signals contaminated by noise, we are 
forced to search for approximate phase synchrony, i.e. 
0i (t) — 02 (t) ~ const. Thus, the studies of synchroniza- 
tion involve not only the determination of instantaneous 
phases of signals but also the introduction of some sta- 
tistical measure of phase locking. 

The wavelet transform is an integral transform for 
which basis functions, known as wavelets, are well lo- 
calized both in time and frequency [20j. Moreover, the 
wavelet basis can be constructed from a single func- 
tion ij;(t) by means of translation and dilation ip(a; t) = 
ip(t—to/ci). The function ip{t) is commonly referred to as 
the mother function or analyzing wavelet. The wavelet 
transform of function s(t) is defined as 



W[s](a,to) = -= / s(t)iP*(a;t )dt, (1) 



where ^*{t) denotes the complex conjugate of ip{t). In 
this work we employ the Morlet wavelet: 



i>{t) = V^T b e 27Tif ^ e - t2/fb 



(2) 



and set the bandwith parameter /& as well as the center 
frequency f c to 1. For a given sampling period 5t, it is 
possible to associate a pseudofrequency with scale a: 
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Obviously, the dual localization of wavelets makes the 
above frequency assignment approximate. 

The instantaneous phase <f>{to) of a signal s(t) can be 
readily extracted by calculating a wavelet transform with 
a complex mother function |lfl.l2J|: 



exp[i(f>(a,t )] = W[s](a,to)/\W[s](a,t )\, (4) 

where we explicitly indicated the dependence of phase 
on the scale a to emphasize that we are investigating 
frequency-specific synchronization, i.e. transient phase- 
locking. 

Following Tass et al. J3| we characterize the strength 
of phase synchronization with the help of the index 7: 



7 = (H ma , - H)/H n 



(5) 



derived from the Shannon entropy H. In the well known 
formula H — J2k=i Pk^ n Pk, N is the number of bins and 
Pk is the relative frequency of finding the phase difference 
within the fe-th bin. Due to normalization in Ipjfl. the syn- 
chronization index lies in the unit interval < 7 < 1. A 
vanishing index 7 = corresponds to uniform distribu- 
tion of phase differences (no synchronization) while 7=1 



corresponds to perfect synchronization (phase locking of 
the two processes). 

We have applied the synchronization theory formal- 
ism (cf. equations {U and ©) to analyze the instanta- 
neous phase difference between ICP and ABP time se- 
ries. Both pressures were averaged over a cardiac cycle. 
The hemodynamic time series were invasively acquired 
during long-term monitoring of patients with traumatic 
brain injuries or spontaneous cerebral hemorrhage. The 
study comprised 10 juvenile patients who were admit- 
ted to the Pediatric Intensive Care Unit of Doernbecher 
Children's Hospital and 10 adult subjects who underwent 
the surgery at the Department of Neurosurgery of Opole 
Regional Medical Center. 

In Table ID we collected the values of mean arterial and 
intracranial pressures for a juvenile patient with trau- 
matic brain injury. From two long-term monitoring ses- 
sions, we have chosen 30 min segments from the time 
series and labeled them A, B and C. 

In the left column of Fig. Q] we present the synchrony 
analysis for data segment A from the first session. Dur- 
ing this period the intracranial pressure remained at an 
elevated but physiologically acceptable level. It is inter- 
esting that for larger scales the phase difference between 
ICP and ABP evolves very slowly. Moreover, the dis- 
tribution of colors indicates that the normalized phase 
difference fluctuates most of the time around 0.5. This 
behavior is in sharp contrast with the phase dynamics for 
segment B from the second session (cf. right column of 
Fig. QJ. During this second period the increase in ICP 
was accompanied by the drop in ABP, which is strong 
evidence of the intermittent failure of regulatory mecha- 
nisms. This failure resulted in insufficient perfusion pres- 
sure. Essentially bichromatic structure of right-top panel 
in Fig. Q] is a clear indication of strong synchronization 
between arterial and intracranial pressure. This observa- 
tion is corroborated by the plot of phase for a = 50 which 
elucidates that the phase variability merely amounts to 
rapid transitions between and 1, equivalent values from 
the point of view of synchronization. 

In Fig. [2] the synchronization parameter 7 for the 
hemodynamic data analyzed in Fig. ^is drawn as a func- 
tion of scale a of the complex Morlet wavelet transform. 
For both curves in this plot, the prominent peak is fol- 
lowed by a plateau. Please note that for the chosen val- 
ues of the Morlet wavelet parameters and St = Is (mean 
cardiac interbeat interval), 1/a gives the approximation 
of the pseudofrequency of the wavelet basis functions 
ip(a;to), cf. ©■ The appearance of such high peaks in 
the low a/high pseudofrequency region is not surprising 
since these maxima are merely the manifestation of the 
inability of cerebral vessels to respond to rapid changes 
in arterial blood pressure. In patients with severely re- 
stricted cerebral blood flow, the characteristic peaks are 
often missing. In Table U we present the values of the 
synchronization strength averaged over scales 30 to 100: 



730:ioo- The maximum value of 0.34 corresponds, as ex- 
pected, to the high pressure/low perfusion episode. To 
assess statistical significance of the observed increase we 
calculated 730:100 for 49 surrogate hemodynamic time se- 
ries |22J. For the surrogate data we found the mean of 
730:ioo equal to 0.07 and since the maximum value was 
0.13, we can reject the null hypothesis that the increase 
was accidental, at the 97.5% level of confidence. 

Fig. shows the instantaneous phase calculations 
for an adult patient after a second massive subarach- 
noid hemorrhage which resulted in severe cerebral edema. 
During monitoring, the average ABP was equal to 81 
Torr and the average ICP was equal to 73 Torr. In this 
case, not only was ICP extremely high but also perfusion 
was negligible. The structure of the phase map (Fig. EJl 
and high value of 7 = 0.30 reflects strong entrainment of 
arterial and intracranial pressure time series. Thus, there 
is no doubt that the phase synchronization observed in 
the two case studies is pathological. In fact, for seven 
patients with a good clinical outcome the average value 
of 7 was 0.09. On the other hand, for five patients with 
severe injury (high ICP, low cerebral perfusion) and poor 
clinical outcome the synchronization parameter was con- 
sistently high and varied between 0.30 and 0.70. 

The relations between the statistical properties of ar- 
terial and intracranial pressure fluctuations are poorly 
understood. Progress in understanding has undoubtedly 
been hindered by the physiological complexity of mech- 
anisms which affect intracranial hemodynamics, see J23J 
for references. In the widely accepted phenomenological 
description of the interaction between ICP and ABP, un- 
der normal conditions a decrease in ABP results in the 
vasolidation of cerebral vessels which increases cerebral 
blood volume and consequently ICP. In pathology, cere- 
bral vessels are non-reactive and changes in arterial blood 
pressure are passively transmitted to ICP. 

Steinmeier et al. and Czosnyka et al. [113, |2J, |25| have 
introduced moving cross-correlation indices which quan- 
tify the reactivity of vessels to changes in ABP. The clini- 
cal studies demonstrated that these indices are positively 
correlated with the high intracranial pressure, low admis- 
sion Glasgow Coma Scale score and poor outcome after 
injury. However, moving correlation or other coherence 
measures based on spectral analysis cannot separate the 
effects of amplitude from those of phase in the interre- 
lations between ABP and ICP signals. Complex wavelet 
analysis provides an effective tool for the investigation of 
phase relations in a chosen frequency range and conse- 
quently can shed new light on inherently nonstationary 
cerebral hemodynamics. 

It is worth pointing out that cerebrovascular vasomotor 
reactivity reflects changes in smooth muscle tone in the 
arterial wall in response to changes in transmural pres- 
sure. The same mechanism underlies cerebral autoregu- 
lation - a fundamental neuroprotective mechanism |ll| . 
In further studies we shall thoroughly test the clinical sig- 



nificance of measures derived from instantaneous phase 
dynamics for assessment of autoregulation integrity. In 
addition, we shall establish the connection between the 
phase synchronization approach introduced in this pa- 
per and those based on moving cross-correlation indices 
[24L |25j or scaling properties of arterial and intracranial 
pressure time series J2(l. l271|. 
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JV 


1 

ABC 


2 
ABC 


ABP 


78.9 77.0 74.6 


59.3 63.3 70.7 


ICP 


13.9 13.0 13.0 


16.4 19.8 9.2 


730:100 


0.13 0.08 0.13 


0.17 0.34 0.13 



Table I: The values of mean arterial ABP and intracranial 
ICP pressures along with the averaged synchronization in- 
dex 730:100 for the 30 min data segments recorded during the 
monitoring of a juvenile patient with traumatic brain injury. 
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Figure 1: Normalized instantaneous phase difference between 
the arterial and intracranial pressure calculated with the Mor- 
let wavelet for 100 integer values of scale a. Left column cor- 
responds to segment A in the first monitoring session, right 
column to segment B in the second (c/. Table |IJ . Bottom 
panels show the time evolution of phases for a — 50. 
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Figure 2: Synchronization parameter 7 as a function of the 
scale a of the complex Morlet wavelet transform. Solid line 
corresponds to segment A in the first monitoring session, 
dashed line to segment B in the second. In Table|l|we present 
the value of 7 averaged over scales 30 to 100 ( as indicated in 
this graph by the horizontal arrow). 
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Figure 3: Normalized instantaneous phase difference between 
the arterial and intracranial pressure for a patient with ex- 
tremely high ICP and negligible perfusion pressure. The cal- 
culations were done in the same way as in Fig. U 



